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We analyze discrete surface modes in semi-infinite binary waveguide arrays, which can support 
simultaneously two types of discrete solitons. We demonstrate that the analysis of linear surface 
states in such arrays provides important information about the existence of nonlinear surface modes 
and their properties. We find numerically the families of both discrete surface solitons and nonlinear 
Tamm (gap) states and study their stability properties. 

PACS numbers: 



Discreteness effects are known to stabilize surface 
modes in nonlinear lattices 0, Q , and stable light self- 
trapping at the edge of a nonlinear self-focusing lattice 
accompanied by the formation of a discrete surface soli- 
ton has recently been demonstrated experimentally [3J 
Such unstaggered discrete surface modes can be treatec 
as discrete solitons trapped at the edge of a waveguide ar- 
ray when the beam power exceeds a certain critical value 
associated with a strong repulsive surface energy U|. 

On the other hand, staggered linear surface modes arc 
known as Tamm states |l5|, and they were first found ir 
solid state physics as localized electronic states at the 
edge of a truncated periodic potential; an optical ana- 
log of linear Tamm states has been demonstrated foi 
the case of an interface separating periodic and homo- 
geneous dielectric media 6*. ^^From that perspective 
staggered surface gap solitons in defocusing semi-infinite 
periodic media, recently introduced theoretically 7] and 
observed experimentally [g , provide a full analogy to lo- 
calized electronic surface Tamm states being an optical 
realization of nonlinear Tamm states. 

The aim of this Letter is twofold. First, we clarify im- 
portant links between different types of linear and non- 
linear surface modes and demonstrate that the analysis 
of linear surface states in waveguide arrays with defects 
provides important information about possible existence 
of nonlinear surface modes. Second, we study discrete 
nonlinear surface states in semi-infinite binary waveg- 
uide arrays, previously introduced theoretically [3, llfl 
and then studied experimentally [I^, and find numer- 
ically the families of both discrete surface solitons and 
nonlinear Tamm (gap) states and analyze their stability. 

We consider propagation and localization of light in 
a semi-infinite periodic binary array of alternating wide 
and narrow weakly coupled optical waveguides, as shown 
schematically in Fig.^top). In infinite binary waveguide 
arrays, the properties of spatial discrete solitons can be 
effectively managed by controlling the geometry of the 
array 9, 10, llj. Following the earlier analysis 9,[l3|, we 
describe the binary array within the tight-binding ap- 
proximation, where the total field is decomposed into a 
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FIG. 1: (a,b) Spectra of the linear surface modes in a binary 
array of alternating thick (A) and thin (B) optical waveguides 
with a surface defect (C) for two configurations of a semi- 
infinite structure shown schematically on top. Shaded areas 
inark two bands of an infinite binary array (p — 0.6). Marked 
points correspond to the examples shown in Figs.2(a-d). 



superposition of weakly overlapping modes of the indi- 
vidual waveguides of two kinds (A and B; wide and nar- 
row). The corresponding equation for the mode ampli- 
tudes takes the form [3| 



.dEr^ 
dz 



PnEn + (En-l + E^+l) + -fn\En\^ E^ = 0, (1) 



where n = 0, 1, . . ., and E_i = due to the structure ter- 
mination. Here p„ characterizes the linear propagation 
constant of the mode guided by the n— th waveguide, 
7„ are the effective nonlinear coefficients (we consider 
Kerr-type medium response) . For the structure shown in 
Fig. da, top), we have p2n+i = -p, P2n = P, and for 
Fig. n^b, top), p2n+i = P, P2n = —p, whcrc p > defines 
the difference between the propagation constants of the 
modes guided by wide and narrow waveguides. 



According to Eq. (Q , the linear Bloch-wave dispersion 
is defined as Kf, — cos~^(— 77/2), where 77 = 2 + p^ — 
/3^. The transmission bands correspond to real Kb, and 
they appear when /3_ < |/3| < /3+, where /3_ = \p\ and 
/3+ = (p^ +4)^/^. The band-gap structure is presented in 
Figs.n^a,b). The upper gap at /3 > /?+ is due to the effect 
of total internal reflection (TIR gap) , whereas additional 
BR and BR2 gaps appear for |/3| < /3_ and (3 < — /9+, 
respectively, due to Bragg reflection. 

In a general case, we assume that the edge waveguide 
(C) differs from both broad (A) and narrow (B) waveg- 
uides, and it is described by a different normalized propa- 
gation constant, po- Then, we study linear surface modes 
for two possible configurations of the binary array with 
the defect at the edge [see top insets in Figs.^a,b)]. We 
look for the localized solutions of Eq. (^ at 7„ = in 
the form £"„ = M„exp(i/3z), where P is the propagation 
constant and u„ is the mode profile. Such eigenmodes 
of the semi-infinite periodic structure with a surface de- 
fect can be found as truncated Bloch waves of the cor- 
responding infinite structure, where the equation for the 
edge waveguide defines an effective boundary condition, 
Po — —pi + [1 + exp(— ii4rb)]/(/3 — pi). The mode be- 
comes exponentially localized inside the gaps where the 
Bloch wave number is complex, and lvaKb{(3) > 0. In- 
deed, we find that spatially localized modes can appear in 
the semi-infinite total internal reflection (TIR) gap [solid 
curve in Figs.n^a,b)] or inside the Bragg reflection gaps 
[dashed and dotted lines in Figs.^a,b) correspond to BR 
and BR2 gaps, respectively]. 

We find that the surface localized modes may ap- 
pear in all three gaps of the wave spectra, as shown 
in Figs. n^a,b), where the marked values of po are: 
p£^ = (p/2) ± ^(p/2)2 + 1, for the structure in Fig. Ufa) 

and pj_^ = —(p/2) ± ^(p/2)2 -I- 1), for the structure in 
Fig. n^b) . Examples of the linear surface modes from dif- 
ferent gaps corresponding to the binary array shown in 
Fig. n^a) are presented in Figs.|2ta-d). 

One of the major observations is that, for the case 
when there is no surface defect and the binary array is 
just cut, i.e. pq — p [see Fig. ^a)] or po = —p [see 
Fig. n^b)] , similar to the original formulation of the prob- 
lem for the Tamm states J5J, no linear localized surface 
modes exist in the system, similar to the case of an array 
composed of identical waveguides discussed earlier 0|. 
The existence of linear surface modes becomes possible 
solely due to the presence of a defect located at the 
edge of the waveguide array, and such localized modes 
are associated with a specific shift of the beam propa- 
gation constant when po 7^ p [Fig. ^a)] and po ^ ~p 
[Fig. n^b)] . The structure of a particular surface mode 
supported by a defect reflects its location with respect 
to the photonic bandgap of the binary array: i.e. it is 
unstaggered in the TIR gap [Fig.j^Ja)], staggered in the 
BR gap [Figs.[2tb,c)], and doubly-staggered in the BR2 
gap [Fig. Hd)]. 

As is well established in many problems of nonlinear 
optics, a nonlinear dependence of the optical refractive 




FIG. 2: (Color online) Examples of the linear surface modes 
in the binary waveguide array shown in Fig. H^a) . The mode 
profiles in (a), (b), (c) and (d) are associated with different 
gaps of the linear spectrum and correspond to the points "a" , 
"b", "c" and "d" marked in Fig. ^a). 



index on the intensity of the incoming light provides an 
important physical mechanism of the intensity-induced 
shift of the propagation constant which can act as an ef- 
fective surface defect. Thus, when we consider the case 
of a nonlinear binary waveguide array, the mode propa- 
gation constant becomes shifted by nonlinearity, and for 
the focusing nonlinear response this shift will be always 
positive. We analyze below two cases where either a wide 
or a narrow waveguide is placed at the edge of the array, 
i.e. when the defect waveguide is created solely by the 
nonlinearity, similar to the case of nonlinear truncated 
lattices studied earlier Q]. Just as in the linear case, we 
look for spatially localized solutions of Eq. ^ now with 
the nonlinear term (we neglect the differences in the ef- 
fective nonlinear coefficients at the A and B sites and 
assume 7„ = 1), in the form £"„ = M„exp(z/?z), where 
/3 is the propagation constant, and the function u„ de- 
scribes the soliton profile. 

In Figs. 01 and 01 we show the families of nonlinear 
surface solitons for two kinds of termination. In the case 
when the edge waveguide is narrow [Fig. (31, the nonlinear 
Tamm states (i.e. nonlinear surface modes in the BR 
gap) appear at low power, in agreement with the linear 
analysis predicting the existence of surface state for a 
small change of the propagation constant [dashed curve 
in Fig. n^b)]. On the other hand, the existence of the 
surface mode in the TIR gap requires the mode power 
to exceed some threshold value, similar to the critical 
value of the defect strength found in the linear theory 
[solid curve in Fig. ^b)]. When the binary lattice is 
truncated at the wide waveguide [Fig. 0] , a finite power 
is required to support the surface solitons in either TIR 
or BR gaps, refiecting the properties of linear TIR modes 
which can be supported only when a positive shift of the 
propagation constant exceeds a certain threshold value 
[solid and dashed curves in Fig.^a)]. 
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FIG. 3: (Color online) (a) Families of nonlinear Tamm (gap) 
states and discrete surface solitons in the binary array trun- 
cated at the narrow waveguide. Solid and dashed curves cor- 
respond to the stable and unstable branches, respectively. 
Schematic of the array is shown on the top. (b,c) Charac- 
teristic examples of the soliton profiles corresponding to the 
points "b" and "c" marked in (a). Detuning parameter of the 
array is p = 0.6. 
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FIG. 4: (Color online) Same as in Fig. |3 but for the binary 
array truncated at the wide waveguide. 



Next, we employ the beam propagation method to 
study the sohton stabihty. In Figs.|3Ja) and^fa), dashed 
branches of the curves indicate unstable surface modes. 
Whereas oscillatory instabilities 1^ may arise for the 
solid branches, we have verified that the solitons corre- 
sponding to the marked points on the sohd curves demon- 
strate stable propagation for more than 100 coupling 
lengths even in the presence of initial perturbations. 



Finally, we mention that the recent studies of discrete 
gap solitons in binary waveguide arrays provided an ex- 
perimental evidence |ll| of a controUed sohton genera- 
tion and power-dependent soliton stee ring in agreement 
with earlier theoretical predictions [^ [l^ . The surface 
modes in the form of discrete surface solitons and nonhn- 
ear Tamm (gap) states described in this paper extend the 
class of such discrete nonlinear modes, and they can be 
treated as the discrete solitons trapped at the edge of a 
waveguide array when the beam power exceeds a certain 
critical value associated with a strong repulsive surface 
energy |j] . We believe that our results open a road for the 
experimental studies of intriguing properties of discrete 
surface gap solitons in engineered periodic structures. 



In conclusion, we have analyzed linear surface states 
at the edge of a semi-infinite binary array of optical 
waveguides, and revealed a link between the linear sur- 
face modes supported by the defects and nonlinear local- 
ized states in the defect-free semi-infinite structure. We 
have found families of both discrete surface solitons and 
nonlinear Tamm states and analyzed their stability. 
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